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ABSTRACT 

The results of infrared photometry from 2 ym-160 urn of AFGL 2591 
12 

and CO observations of' its associated molecular cloud and high 
velocity molecular outflow are presented and discussed. Observations 
between 60 ym and 160 ym provide the first detailed far-infrared 
spectrum of this source. The observed luminosity between 2 and 160 ym 

4 

is 6.7 x 10 L at a distance of 2 kpc. We estimate the total 

® 

4 

luminosity to be ^ 9 x 10 L q . The 2ym to 160 y.m spectrum of AFGL 

2591 is interpreted in the context of a model in which a single embedded 

object is the dominant source of the infrared luminosity. We determine 

this object to be surrounded by a compact, optically thick dust shell 

with a temperature in excess of several hundred degrees kelvin. We 

estimate the extinction to this source to be between 26 and 50 visual 

magnitudes. Our data is consistent with earlier suggestions that the 

object which powers the near-infrared source also is responsible for 

ionizing the compact HII regions near AFGL 2591. We have determined the 

absolute position of the infrared source at 10 ym to an accuracy of ± 

1" . This indicates for the first time that the IR source and H^O source 

are not coincident. 

12 

Our CO observations show the high-velocity molecular flow near 

AFGL 2591 to be extended, bipolar and roughly centered on the infrared 

emission. The observations suggest that the red-shifted flow component 

extends beyond the boundary of the ambient cloud within which AFGL 2591 
12 

is embedded. The CO observations also show that AFGL 2591 is embedded 
in a molecular cloud with an LSR velocity of -5 km s"^. Comparing this 
velocity with the velocity of H 2 0 maser emission -22 km s" 1 ) suggests 
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that the maser source is part of the blue-shifted high-velocity 
molecular outflow. Its close proximity to the infrared source suggests 
that the origin of the molecular outflow may be very close to (within 
1") the surface of AFGL 2591. 



I. INTRODUCTION 


The observation that most stars, in their earliest stages of 
stellar evolution, produce cold energetic outflows of molecular gas 
(e.g., Bally and Lada 1983) suggests a revision of our understanding of 
the nature of young embedded infrared sources (e.g., Wynn-Wil 1 iams 

1982) . It is now clear that many objects previously thought to be 

quietly accreting protostars are actually the centers of very energetic 
mass explusion. One such object, AFGL 2591, typical of the 

proto-stellar sources, has recently been found to be associated with a 
high-velocity molecular outflow (Bally and Lada 1983, Torrelles et al 

1983) . AFGL 2591 is a very compact (% 0.2", Howell, McCarthy and Low 

4 

1981) and luminous (L £ 3 x 10 L > Merill and Soifer 1974) near- infrared 
source. It is closely associated with an H^O maser source 
(Wynn-Wil 1 iams 1977 ) and also near at least two, weak compact HII 
regions (e.g. Simon et al 1981). Although very luminous and compact, 
there is no evidence for ionized emission at the location of the 
near-infrared source (Thompson and Tokunaga 1979, Simon et al 1981), 
suggesting that the source has not yet evolved to the stage of emitting 
significant ionizing flux. Despite extensive near-infrared studies, 
little is known about the far-infrared spectrum of this object and its 
total luminosity has not been accurately determined. The recent 
detection of high-velocity molecular flow around this object and the 
lack of a good far- infrared spectrum of this object motivated us to 
obtain further observations of this source and reconsider its nature. 

In particular we were interested in determining the characteri sties of 
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the molecular flow and its relation to the infrared source, compact H I I 
regions and water maser. We also were interested in obtaining a good 
far-infrared spectrum in order to obtain a reliable total luminosity and 
investigate the detailed energetics of the source. 

To accomplish these objectives we obtained new mill imeter-wave CO 
and infrared observations of AF6L 2591. The results of these 
observations, which include the detection of a bipol ar molecular flow 
and the first accurate far-infrared spectrum of this source, are 
presented in this paper. We also present new 2-20 pm broadband 
spectroscopy of the near-infrared source and a 10 pa absolute position 
measurement of the source, which significantly improves upon previous 
determinations made at shorter wavelengths. The telescopes and 
observing techniques employed for this study are described in Section 
II, the results and analysis of observations in Section III and the 
implications of the results in Section IV and Section V. 
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II. TELESCOPES AND OBSERVING TECHNIQUES 
a ) Millimeter-wave Spectroscopy 

1 o I r 

Observations of the J=l-0 transition of C O at 2.6 mm were 

initially obtained in 1981 February with the 11-meter mill imeter-wave 

2 

telescope of the National Radio Astronomy Observatory (NRAO) in 


? 

The National Radio Astronomy Observatory is operated by the Associated 
Universities, Inc., under contract with the National Science Foundation. 


Arizona. Subsequently the majority of the observations reported here 
were obtained in 1983 June with the new 12-meter surface of the NRAO 
millimeter-wave telescope. Both sets of CO observations were made with 
the same cooled, dual-channel, mixer receiver which was operated in 
conjunction with a 512-channel filter bank spectrometer, which provided 
spectral resolutions of 0.25 MHz and 0.50 MHz. The filters were split 
so that each polarization channel of the receiver was coupled to two 
independent 128-channel banks of 0.50 MHz and 0.25 MHz filters 
respectively. Data taken at the same resolution in the two receiver 
channels were averaged together. All spectra were observed in the 
position switching mode using the same "off" position (i.e. 
a = 20 h 2 m 35 S .9, 6 = 43° 01' 16"; 1950.0) which was pre-determined to be 

free of significant CO emission. Data taken with both telescopes were 
initially calibrated using a rotating chopper wheel following U1 ich and 
Haas (1976). The 11-m data were finally calibrated by observing Orion 
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•k 

and IRC +10° 216 and adopting temperature scales such that = 60.0 K 

k 

and = 4.1 K for these two objects respectively. Final calibrations 
for the 12-m data was accomplished by observing M17 SW and IRC +10° 216 

ic : k 

and adopting the observed temperatures = 35.2 K and = 5.8 K for 

these two objects respectively. The M17SW temperature agrees within 10% 

* 

with standard 11-m value of = 39.0 K (Thronson and Lada 1983, Ulich 
and Haas 1976) while the value for IRC +10° 216 is 38% higher than the 
standard 11-m calibration. This is an expected result of the higher 
aperture efficiency of the 12-m surface. 

b) Near- and Mid- Infrared Photometry 

AFGL 2591 was observed during the autumn of 1982 using a Ge(Ga) 
bolometer system on the 2.3 m telescope of the Wyoming Infrared 
Observatory (WIR0). Observations were obtained with broadband filters 
centered at the wavelengths listed in Table 2 and used a 6.6" aperture 
and a chopping secondary with a 30" north- south throw. The filter 
characteristics and absolute calibrations are the same as presented by 
Gehrz, Hackwell and Jones (1974). The internal uncertainty for these 
observations was less than 2% in flux (1^ RMS). Therefore, the total 
uncertainty is dominated by systematic errors in the absolute 
calibration of standard sources. We estimate this uncertainty to be ±5% 
at 2.3 ym, 3.6 ym and 4.9 ym, and ±10% at the longer wavelengths. 



8 


e) Far- Infrared Photometry 

AFGL 2591 was observed during 1982 August at far-infrared 
wavelengths using the 0.9 m telescope onboard NASA's Kuiper Airborne 
Observatory (KAO). A mul ti-channel , multi-filter Ge(Ga) system was used 
to map the source over an area of 3' x 6' (east-west by north- south) . 

II 

The angular resolution of the system was 49 . Four filters with 

effective center wavelengths of 60 pm, 95 ym, IlOym, and 160 ym were 
used to obtain far-infrared flux densities. The absolute calibration 
uncertainty in the flux density is approximately ±20%. The statistical 
uncertainty was about ±1% (la RMS) and therefore can be ignored. 
Details concerning filter characteristics and calibration procedures 
have been described elsewhere (Thronson and Harper 1979, Loewenstein 
et al . 1977). The W3(0H) source was used for calibration. Data were 
taken by chopping the secondary to a reference position 7.5' north and 
south of AFGL 2591. The maximum far-infrared emission was located at 
the 10 ym peak, within the ±15" positional uncertainty, and the object 
was unresolved. Weak, extended emission was observed from the source 
throughout the region mapped, but at a level 0.5-1% of the peak channel. 
It is unlikely that this flux contributes more than n, 20% to the total 
measured at the peak of AFGL 2591. 
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III. RESULTS AND ANALYSIS 

a ) CO Observations 

Spectra of D 0 emission were obtained at 21 locations in a grid 

centered on AFGL 2591. Ten positions were observed with the 11-m 

telescope and 14 positions with the 12-m telescope, while three 

positions were observed in common with both telescopes. At the infrared 

peak of AFGL 2591, the spectrum consists of two strong but narrow 

emission line components superposed upon broad emission wings which 

extend over about 40 km s' 1 (e.g. Bally and Lada 1983, Torre! 1 es et al 

1983). The strongest narrow line component is at a velocity of about -5 

km s' 1 , and the weaker component at about +1 km s' 1 . The temperatures 

* -1 

of the two components centered on AFGL 2591 are T^ = 15.2 K (-5 km s ) 

★ — I 

and T^ = 5.1 K (+1 km s ) when observed with the 11 m telescope and 
T^ = 18.5 K (-5 km s"' 1 ') and T^ = 5.4 K (+1 km s' 1 ) when observed with 
the 12 m telescope. The difference in the two sets of observations is 
apparently a result of the fact the +1 km s' 1 component is considerably 
more extended than the -5 km s 1 component. The higher aperture 
efficiency of the 12 m surface will result in a 1 arger antenna 
temperature for less extended objects. In fact, at the positions 
observed in common the antenna temperatures of the +1 km s' 1 component 
agree to within 10% while the -5 km s" 1 component is systematically 
stronger by about 20% in two positions (0,0 and 1W,1S). This suggest 
that emission from -5 km s" 1 does not fill the beam uniformly and is 
clumpy. To make a map of the distribution of emission of the -5 km s' 1 
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component, we scaled the 12-m data to match the 11-m data at the 
position of AFGL 2591. The results are shown in Figure 1 in the form of 

* _i 

a contour map of T^ (-5 km s ). The positions where spectra were 

observed are also indicated on the map by crosses. We find that the 
12 -1 

CO emission (at -5 km s ) peaks on AFGL 2591, which is also apparent 
when both sets of data are treated separately. The clear intensity 
enhancement of the -5 km s"* cloud at the position of AFGL 2591 suggests 
that the infrared source is embedded in this cloud. Emission at -5 km 
s" 1 falls off rapidly to the northeast of AFGL 2591, but more gradually 
to the southwest. On the other hand, there appears to be little 
variation in the intensity of the 1 km s" 1 component across the mapped 
grid. This can be seen in Table 1 where we list the observed antenna 
temperatures of the 1 km s" * component at each position of the grid. At 

most locations the antenna temperature of this component is within 0.5 K 

* -1 
of T^ = 5.0 K. This confirms our earlier suggestion that the 1 km s 

cloud is extended and uniform. 

In order to investigate the distribution of high-velocity CO 
emission around AFGL 2591, we determined the integrated intensities in 
the red and blue high-velocity wings for CO profiles at each position in 
the grid as follows: 


RED 


14.5 

* 

V* 

3.8 


and 


BLUE = 


- 10.2 

* 

* 

t a dv 

• -21.9 
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These integrated intensities are also listed in Table 1. The maximum 

values occur near AFGL 2591 and in the blue component, which at the 

infrared source is four times stronger than the red component. However, 

there is considerable variation of the ratio of BLUE/RED across the 

source and our observations clearly show the flow in AFGL 2591 to be 

bipolar. The approximate boundaries of the regions where the integrated 

intensities in the red and blue components exceed 3 K-km s"' 1 ' are 

indicated in Figure 1. Again the bipolar nature of this flow is 

apparent. Figure 2 shows spectra obtained at two positions in the flow 

marked with x's in Figure 1. They dramatically illustrate the distinct 

spatial separation of the red and blue high-velocity flow components. 

Indeed, the most intense high-velocity emission we measured in the red 

flow occurs at the edge of our grid, apparently well beyond the boundary 

of the -5 km s _1 cloud. Unfortunately we did not observe enough 

positions to determine the true extent and distribution of the 

13 

high-velocity gas. In addition, we did not obtain CO observations and 
we are unable to reliably estimate either the masses or the energetics 
for the high-velocity flow. Recent observations by Torrelles et al 
(1983) of a similar sized region around AFGL 2591 with the 11-m 
telescope did not indicate a bipolar nature to the high-velocity 
emission. However, their data did not have sufficient signal -to-noise 
to detect the red component of the flow and were only sensitive enough 
to map the distribution of the brightest high-velocity blue-shifted gas. 
Their map of the blue component of the high-velocity flow is in 
satisfactory agreement with our results. 
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b) Infrared Observations 


The results of our infrared photometry from 2 to 160 ym are 
presented in Table 2. The resulting spectrum is plotted in Figure 3. 
Our 2-20 ym spectrum agrees well with previously published observations 
(Merrill and Soifer 1974; Wi liner et al 1982). The spectrum is very 
similar to "protostars" such as BN, W3 IRS 5 and NGC 7538 IRS 9, objects 
all associated with high-velocity molecular flows. Integrating under 
the observed spectrum we obtain the first accurate luminosity for this 
source. We find L/L = 6.7 x 10 (D/2 kpc) , where D is the distance 

to AFGL 2591 in kiloparsecs. Roughly half the luminosity is observed 
between 2 and 20 ym, as the figure shows. 

The 2-20 ym spectrum of AFGL 2591 mimics very closely the shape of 
a blackbody, excluding the 10 ym "silicate" feature. By fitting a 
Planck distribution to the observed points, we can obtain estimates of 
the temperature of the warm dust and the angular size of the emitting 
region. We find that the observed points fit very closely a blackbody 
distribution given by a temperature, T^ = 460 K, and a source 
diameter, 0 = 0."13. This source diameter is consistent with measured 
upper limits to the source sizes obtained by speckle 
interferometry 0 < 0."2 at 2.2 ym; (Howell, McCarthy, and Low 
1981); 0 < 0."5 at 10 ym (D. Mozurkewich, private communication). 
However, the derived dust temperature and source size are actually lower 
limits to the true values since the 10 ym absorption feature suggests 
that the infrared source is heavily extinguished. 
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The far- infrared spectrum can likewise be fit fairly well by a 
blackbody spectrum, but this is due in a large part to the small number 
of filters we used. Figure 3 shows a best fit, assuming that the 
intrinsic spectrum is, in fact, a blackbody. The resultant values (T Q = 
73 K, with 0=8") are actually lower limits to the true values. For 
example, most far-infrared spectra are assumed to arise from 
optical 1 y- thin isothermal dust emission and that the source spectrum is 

of the form F « v B (T). If this is the case for AFGL 2591, then T n ^ 

v v U 

110 K. However, Thronson and Harper (1984) have emphasized that 
overlying, non-emitting dust can alter a far-infrared spectrum and make 
it appear significantly "cooler" than the true dust temperature. In 
addition to this, Thronson and Harper further argue that for the vast 
majority of sources observed to date, there is not sufficient spectral 
coverage to determine with certainty whether or not this effect is 
taking place. Since this is the case for our data on AFGL 2591, we are 
limited to saying that the dust that dominates the far-infrared emission 
is characterized by T Q > 73 K, although more likely, T Q ^ 100 K. 

The blackbody fits to the spectrum of AFGL 2591 enable us to 
estimate the total luminosity of the source including the contribution 
from enission outside the observed wavelength range. Integrating under 

4 

the two black body curves, we estimate a total luminosity of 9 x 10 
2 

(D/2 kpc) L 0 for the source. 

Finally, the excellent pointing characteristics of the WIRO 2.3 m 
telescope enabled us to measure the absolute position of the N band 
(10 pm) emission to an accuracy of ±1" (Table 2). This represents a 
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substantial improvement over previous position estimates of AFGL 2591 
made at 3.5 ym with ±2" accuracy ( Wynn- Will iams et al . 1977). The 
location of the 10 ym source is shown in Figure 1 relative to the 
locations of the H 2 0 maser associated with AFGL 2591 and continuum 
emission from compact HII regions in the vicinity (Wynn-Will iams et al . , 
Simon et al . 1981). Our improved position determination shows that the 
maser spot and IR source are not coincident, as previously thought. The 
infrared source is significantly displaced from the compact HII regions 
as first determined by Wynn-Will iams et al . (1977). We note that our 
simple interpretation assumes that 3.5 ym emission arises from the same 
material as does the 10 ym emission. 
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IV. DISCUSSION 


a) Nature of AFGL 2591 

As shown in Figure 1 the immediate vicinity of AFGL 2591 contains 
two compact radio continuum sources and an H^O maser spot, as well as 
the compact near-infrared source. In addition, two smaller, less bright 
compact continuum sources are also known to be contained within this 
region (Campbell 1984). The compact continuum source southwest of the 
infrared source has been studied by many authors (Wendker and Baars 
1974, Wynn- Will jams et al . 1977, Simon et al . 1981, Brown 1974, Campbell 
1984) and although the fluxes determined between 2-15 GHz are not in 
particularly good agreement, they are consistent with a spectrum of an 
optically thin, photoionized HII region. If all four continuum sources 
are optically thin H II regions within the cloud and are photoionized by 

4 

individual ZAMS stars they would contribute about 2 x 10 L @ to the 
total far-infrared luminosity in the region. That is about 25% of the 
total luminosity we estimated for the region from our far-infrared 
observations. However there appears to be no strong 2 yn or 10 ym 
sources at any of these positions, and it is likely that these compact 
continuum sources are not powered by stars which contribute 
significantly to the far-infrared luminosity we observed. Indeed, 
Thompson and Tokunaga (1979) have suggested that the ionization of the H 
II region southwest of the infrared source is provided by UV radiation 
escaping from a hole in an otherwise thick dust shell around AFGL 2591. 
Such a model also provides a natural explanation for the presence of the 
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reflection nebula located approximately 20" southwest of AFGL 2591 on a 

line connecting the H II region and infrared source (Kleinmann and 

Lebofsky 1975). To provide the requisite amount of photoionization for 

the compact H II region, AFGL 2591 must produce an ionizing flux 

equivalent to that of at least an 07.5 ZAMS star. Since the luminosity 

4 , 

of an 07.5 star is about 8 x 10 L g (Panagia 1974), our estimated total 

4 

luminosity for AFGL 2591 of 9 x 10 L q is consistent with the idea that 
AFGL 2591 could be the source of inoization for at least the largest 
compact continuum source in the region. Our observations would be 
consistent with this model as long as the distance to AFGL 2591 is 
greater than 1 kpc. Clearly, our observations support the suggestion 
that AFGL 2591 is the dominant energy source in the region. 

As the dominate energy source, AFGL 2591 would also be the driving 
source of the bipolar molecular flow. Since it is believed that such 
outflows are driven by winds emmanating very close to the surface of the 
central young stellar objects (e.g., Lada 1984), it is not surprising 
that holes or channels exist in a dust shell which would permit UV 
radiation to escape and ionize the compact HI I region as well as produce 
the observed reflection nebula. Indeed, reflection nebula are found 
associated with many mol ecu! ar outflows (Strom 1983). The existence of 
weaker H 1 1 regions near AFGL 2591 (see Figure 1, Simon et al . 1981 and 
Campbell 1984) may indicate that the dust shell is beginning to break 
up. However, it is also possible that the radio continuum emission is 
being produced by collisional ionization in a stellar wind (Krolik and 
Smith 1981). This certainly appears to be the case for a number of 
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compact HII regions associated with other high-velocity outflows (e.g., 
Simon et al . 1981; Cohen, Bieging and Schwartz 1982; Bally and Predmore 
1983; and others). Although the existing radio spectra of AFGL 2591 
are not presently suggestive of such a situation, high spectral 
resolution infrared recombination line observations might be able to 
distinguish between an outflowing col 1 is ional 1 y ionized wind and a 
photo-ionized HII region. In addition sensitive searches for 
near-infrared sources at the site of the strongest continuum source 
should easily uncover a BO. 5 star if one was actually ionizing the HII 
region. At 2.2 ym a BO. 5 star at the distance of AFGL 2591 would have 
an apparent magnitude of +10 if extinguished by 50 visual magnitudes. 
However if such a star was present we should have detected it at 10 u m. 

If we assume that the compact near-infrared object is the source of 
all the observed far-infrared emission, we can place interesting 
constraints on its nature. As pointed out earlier, the presence of a 
deep 10 ym silicate absorption feature in the spectrum indicates that 
between 2-20 ym we are observing heavily reddened emission from hot 
dust surrounding a recently formed star. An upper limit to the 
extinction to this hot dust can be estimated by de-reddening the near 
infrared observations and then requiring that the resulting total 

4 

luminosity between 2-20 ym not exceed 9 x 10 L q . If we assume a 
standard interstellar reddening law (e.g., Van de Hulst curve no. 15 
[Johnson 1968]) we find an upper limit to the extinction of about 50 
visual magnitudes between 2.3 ym and 4.9 ym. We can obtain a lower 
limit to the extinction from the depth of the 10 ym absorption feature. 
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Rieke and Lebofsky (1984) have found that = 16.6 ± 2.1 x T . from 
careful measurements of a number of heavily obscured sources. Merrill 
and Soifer (1974) estimated an optical depth of = 1.8 for the 
silicate feature in AFGL 2591. This corresponds to A y = 29.9 ± 3.8 
magnitudes. This is a lower limit to the total extinction, since there 
is likely to be emission as well as absorption present at 10 pm. 
Therefore the observed depth of the silicate feature underestimates the 
actual extinction to the source. We conclude that AFGL 2591 is obscured 
by 26-50 magnitudes of visual extinction. 

The presence of such a large amount of extinction towards the 
near-infrared emitting dust also suggests that the circumstel 1 ar dust is 
quite hot with tempertures in excess of 500 K. The compact nature of 
the near-infrared source further suggests that the range in temperature 
of the hot circumstell ar dust may be small. For example, let us assume 
that the circumstellar dust has an intrinsic energy distribution of a 
single temperature blackbody between 2 and 5 pm. If we require the 
total luminosity of this blackbody to equal that of AFGL 2591, then, in 
principle, we can uniquely determine its temperature and size, once we 
specify the flux density of the blackbody distribution at any single 
wavelength. Because all the observed emission between 2 and 5 pm 
suffers considerable extinction we are uncertain of the intrinsic flux 
at any given wavelength. However we know that the extinction to the hot 
dust is between 26 and 50 visual magnitudes. Therefore if we assume an 
extinction law and a single value of visual extinction we can de-redden 
the 2-5 y m observations and attempt to fit the de-reddened points by a 
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blackbody whose total luminosity equals that of AFGL 2591. For possible 
extinctions between 26 and 50 visual magnitudes and an extinction law 
given by curve #15 of Van de Hulst (Johnson 1968), we find that the 
de- reddened observations can be fit by a single temperature blackbody 
of appropriate luminosity only when A y £ 43 magnitudes. The resulting 
blackbody curve is characterized by a temperature of about 800 K and a 
size of about 0.06 arc seconds. Figure 3 shows the curve fit to the 
de- reddened observations between 2 and 5 urn. 

It is of interest to compare the blackbody curve with our 8-13 pm 
observations. Using the extinction law derived by Rieke and Lebofsky 
(1984) between 8 pm and 13 ym we can predict the fluxes in this 
wavelength range expected from silicate absorption of the underlying 
blackbody emission if A y £ 43 magnitudes. We find that the predicted 
fluxes in the silicate absorption feature are about 1.5 times smaller 
than that actually observed between 8.0 ym and 11.4 ym. In other 
words it appears that the extinction at 10 ym is equivalent to only 
about 23 visual magnitudes. However this may still be consistent with 
our single temperature dust shell model, since we expect the silicate 
feature to have a significant emission component. Longward of 13 ym 
the observed fluxes are greater than that of the blackbody curve and at 
these wavelengths our assumption of dust emitting at a single 
temperature breaks down. In any event, our observations are consistent 
with and perhaps best represented by a simple model for AFGL 2591 in 
which a single object surrounded by a hot (T^ £ 800 K) , compact 

II 

(0 ^ .06 ) dust shell is the source for all the observed far-infrared 
emission . 
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b) The Nature of the Molecular Outflow 


Our CO observations are not comprehensive enough to enable as 
thorough an analysis of the physical properties of this molecular 
outflow as has been possible for a number of other well studied sources 
(e.g., NGC 2071, Bally 1982; AFGL 490, Lada and Harvey 1981; AFGL 961, 
Lada and Gautier 1982). However, despite these limitations our 
observations do show AFGL 2591 to be associated with one of the more 
interesting of the known examples of high-velocity molecular outflow. 

We will now consider briefly some of the more interesting aspects of the 
outflow revealed by our data. 

First, we find the flow to be bipolar, and as indicated by Figure 
2, the spatial separation of red and blue high-velocity components is as 
distinct as has been found for any other flow observed to date. 
However, our observations are not sufficiently extensive to determine 
the full extent or degree of coll imation of the flow. Clearly further 
mapping is needed to address these particular issues. 

Second, we find that the vel oci ty extent of the high-velocity gas 
is as great in the outermost regions of the flow as it is at the 
position of the infrared source at the apparent center of the bipolar 
flow. Thus, there is no evidence for any flow deceleration with 
distance from the central source. In fact, the velocity extent of the 
high-velocity emission is somewhat higher in the outer regions than near 
the center of the grid. Such kinematical behavior is not so clearly 
evident in any other extended source we have studied. Acceleration, or 
lack of deceleration, of high-velocity gas would be expected if the 
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force accelerating the flow was more or less constantly applied over the 

dynamical time scale and if the flow was propagating into a medium with 

a decreasing density gradient. Of particular interest in this regard is 

the extent of the red-shifted flow beyond the apparent boundary of the 

-5 km s ambient cloud as seen in Figure 1 and discussed in Section 

IIA. Although, as can be seen Figure 2, there js_ weak emission 
* _1 

(T^ ^ 1.4 K) at -5 km s in the outer red-shifted region it is not 

clear whether or not this emission in fact originates from ambient gas 

in the outermost regions of the -5 km s - ^ cloud. The intriguing 

possibility remains that in these more extended regions even emission 

observed at -5 km s"* arises in the molecular outflow rather than from 

13 

ambient material. Observations of CO emission are necessary to 

13 

investigate this possibility. Even if CO observations were to show 
that emission at -5 km s~^ arose in the quiescent part of the cloud, it 
is still clear from Figure 1 that the red-shifted flow component has 
broken out into the outermost, low density regions of the molecular 
cloud. It is possible that this flow component is freely expanding and 
that its extent is considerably larger than suggested by our map. 
Recently, Schwartz, Smith and Waak (1983), have suggested that the 
red-shifted flow associated with NGC 1333 outflow is expanding into a 
very low density region of that cloud. However the AFGL 2591 flow may 
be the first molecular outflow which appears to have expanded beyond the 
boundaries of its associated molecular cloud. 

Third, we find the peak integrated intensity in the blue flow to be 
considerably greater than that in the red flow. This is similar to the 
situation for the S 140 outflow (Lada and Wolf 1983) but not typical of 
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outflow sources which generally have their strongest emission in the red 
component (Bally and Lada 1933). Although it is not all clear what 
produces such asymmetries in the flow line profiles, the apparent 
enhancement of blue emission in AFGL 2591 could be due to the fact that 
the blue flow is propagating into a higher density region than the red 
flow, sweeping up more material in a smaller area and appearing 
brighter. In this case we would expect the blue flow to be considerably 
more confined than the red flow. Future mapping could reveal if this 
was the case. 

Fourth, our data enables us to make crude estimates of the flow 
extent and dynamical age. From Figure 1 we estimate the flow to be at 
least 5‘ in extent, which corresponds to a size of 3 pc at a distance of 
2 kpc. This, coupled with the observed velocity extent of the flow 
('u 40 kms -1 ), suggests a dynamical time scale of ^7 x 10^ years (at 2 
kpc distance). These parameters indicate that the AFGL 2591 flow is one 
of the most extended and perhaps evolved of the known molecular outflow 
objects. AFGL 2591 may also be the most luminous driving source of the 
known molecular flow sources. It is perhaps surprising that such a 
luminous young stellar object has not had time enough to evolve to the 
point where it has produced a significant compact HII regions, yet has 
had enough time to generate one of the most extended outflows known. 
Consider that the molecular flow associated with IRc 2 in Orion (the 
only other well studied outflow source of comparable luminosity to AFGL 
2591) has an extent of 0.2 pc and a dynamical age of MO years. 

The lack of a developed HII region around AFGL 2591 might be due to 
the fact that the exciting star is surrounded by an extraordinarily 



dense and dusty shell which has severely retarded the expansion of an 
HI I region. However, the presence of weak continuum sources near AFGL 
2591 suggests that the dust shell is breaking apart as discussed 
earl ier . 

It is also possible that AFGL 2591 is closer than 2 kpc. At a 

4 

distance of 1 kpc it would have a luminosity of 2 x 10 L , a flow 

4 

extent of ^1.5 pc and flow dynamical age of 3 x 10 years, typical of 
the majority of luminous molecular outflow sources. 

Finally, we comment on the relation of the H^O maser to AFGL 2591 
and the molecular flow. Our data indicate that AFGL 2591 is embedded in 
the -5 km s"^ cloud. Moreover it is not coincident with the H^O maser 
source as previously thought. The velocity of maser emission of -21 km 
s’" 1 and -23 km s" 1 (Wynn-Will iams et. al 1977) clearly shows the maser 
source within the blue-shifted high-velocity flow. It is apparently a 
dense clump of gas accelerated to high-velocity by the outflow. Its 
close proximity to AFGL 2591 again suggests that the origin of the 
outflow occurs very near the infrared star. 
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V. SUMMARY 

The primary results of our infrared and millimeter-wave study of 
AFGL 2591 can be summarized as follows: 

1) The luminosity observed between 2 m and 160 m for the AFGL 

4 2 

2591 source is found to be 6.7 x 10 (D/2 kpc) L , suggesting a total 

0 

4 

luminosity at all wavelengths of about 9 x 10 l for a source distance 

0 

of 2 kpc. 

2) The near- and far-infrared observations are consistent with and 
best explained by a model in which a single embedded object is the sole 
source of the observed near- and far-infrared luminosity. This object 
appears to be surrounded by a compact, optically thick dust shell with a 
temperature of % 800 K and an angular diameter of ^0.06", although a 
significant range in temperatures probably actually exists. The 
extinction to this source is estimated to be about 43 visual magnitudes. 

3) In addition, our data are consistent with the suggestion of 
Thompson and Tokunaga (1979) that AFGL 2591 is the source of ionization 
for the compact HI I regions in its vicinity and the source of scattered 
light for a nearby reflection nebula. As the dominant energy source in 
the region, AFGL 2591 is probably the sole source driving the high 
velocity molecular outflow associated with it. 

4) Our CO observations indicate that the high-vel ocity mol ecul ar 
outflow associated with AFGL 2591 is extended and bipolar in nature. 

5) The red-shifted component of the bipolar outflow appears to 
extend beyond the boundary of the ambient cloud in which AFGL 2591 is 
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embedded. The lack of significant deceleration of the red-shifted flow 
may be consistent with its propagation into a regime of low cloud 
density. The blue-shifted component of the flow is considerably 
brighter than the red-shifted flow in the vicinity of the infrared 
source, and appears to be within the cloud boundary. However, it also 
shows no significant deceleration with distance from the flow 

center. The material giving rise to both velocity components may be 
undergoing constant acceleration. 

6) CO observations also indicate that the molecular cloud within 
which AFGL 2591 is embedded has an LSR velocity of -5 km s" 1 . The 
velocities (i.e., -21 and -23 km s"^) of the H^O maser source associated 
with AFGL 2591, suggest therefore that the maser source is actually part 
of the blue-shifted high-velocity molecular flow. 

7) The position of the near-infrared source has been determined to 
an accuracy of ±1" at 10 urn and we find that the infrared source is not 
coincident with the H^O maser spot as previously thought. However, the 
close proximity of the maser spot to the infrared source does suggest 
that the high-velocity molecular flow originates very close to the 
embedded star. 
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AFGL 2591 Source Parameters 


a(1950) 

6(1950) 


20 h 27 m 35^3 ± 0^1 
40° 1' 14" ± 1" 


Infrared Photometry 


WIRO (6"6 Aperture) 9 


2.3 

Pm 
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70 
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Luminosity (2-200 pm, 2 kpc) 6.7 x 10^ L q 
T d (40-200 pm) = 73K; 0 = 8" 

Tj-j (2-20 pm) = 460K (observed); 0 = 0.13 ’ 

Tp (2-20 pm) = 800K (de-reddened) C ; 0 = 0.06 


Ay = 43 mag c . 


The statistical uncertainty of the flux densities is in all cases bette 
than 2%. The calibration uncertainty is estimated to be ±5% at 2.3 pm, 
3.6 pm, and 4.9 pm. At the longer wavelengths, we estimate a ±10% 
calibration uncertainty. 

he estimate a total uncertainty of ±20%. 
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FIGURE CAPTIONS 

Figure 1 

*12 -1 
Contour map of peak ( CO) emission from the -5 km s molecular 

cloud near AFGL 2591. Contour intervals are in steps of 2 K starting 

with the outermost contour which is 4 K. Also plotted are the 

approximate boundaries of the red- and blue- shifted components of the 

high-velocity molecular flow, showing the apparent bipolar nature of the 

12 

flow. The positions where CO observations were taken are indicated by 
"+" symbols. The location of the most intense integrated emission in 
the red-shifted flow is indicated by "+". The "x" symbols mark the 
location where the spectra shown in Figure 2 were obtained. The insert 
shows the small-scale structure of AFGL 2591. The 10 yin peak is 
presented in relation to the H^O maser position (Wynn-Will jams et al . 
1977) and radio continuum emission (adapted from Simon et al . 1931). 

Figure 2 

12 

The CO spectra observed at the two locations indicated in Figure 
1 by the "x" symbols. The distinct assymetry of the red and 
blue-shifted high-velocity line wings is clearly evident, indicative of 
a bipolar flow pattern. Also evident are the quiescent ambient clouds 
at V, SR £ +1 km s' 1 and -5 km s" 1 . 

Figure 3 

The observed infrared spectrum of AFGL 2591, along with simple 
models of the continuum emission at these wavelengths. The observed 
points are connected by straight lines. Our suggested models for the 
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mid-infrared and the far-infrared emission are the two Planck functions 
labeled by temperature and angular diameter. Three shorter wavelength 
points are shown after correction for A y = 43 mags. The total 
uncertainties in each flux density measurement are shown as ±la error 
bars, if they are larger than the symbol. 
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